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Abstract 

Discrete VT — symmetric square wells are studied. Their wave functions are found 
proportional to classical Tshebyshev polynomials of complex argument. The compact 
secular equations for energies are derived giving the real spectra in certain intervals 
of non-Hermiticity strengths Z. It is amusing to notice that although the known 
square well re-emerges in the usual continuum limit, a twice as rich, upside-down 
symmetric spectrum is exhibited by all its present discretized predecessors. 
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1 Introduction 



Undoubtedly, solvable models represent one of the most important and inspiring 
sources of insight in the physical properties of real quantum systems [1]. For illus- 
tration, one may recollect how the concept of the so called shape-invariant analytic 
potentials V(x) [2] opened the path towards a better understanding of supersymme- 
try [3]. Recently, a similar major progress in our understanding of some less standard 
structures admitted by the standard Quantum Mechanics [4] has been initiated by a 
few numerical as well as non-numerical analyses [5, 6, 7] of certain complex potentials 
exhibiting another interesting invariance property, 

V(x) = [V(-x)T • (1) 

It has been argued [7] that the latter invariance may prove highly relevant in physics 
as it makes all the Hamiltonian VT— symmetric, i.e., commuting with the product of 
parity V and complex conjugation T while, in its turn, the latter operator represents 
the time-reversal operation in the language of phenomenological considerations. 

In the resulting slightly innovated formalism of quantization (called, for defmite- 
ness, VT — symmetric Quantum Mechanics, PTSQM [8, 9]), you admit that Hamil- 
tonians need not necessarily be Hermitian in the usual trivial metric I in Hilbert 
space, H ^ W . As an illustrative example one may recollect the ordinary differen- 
tial Schrodinger equation on a finite interval, 

^(x)=E^(x), V(±1) = (2) 

where PTSQM admits potentials V(x) of eq. (1) which are not real but merely 
"parity plus time-reversal" - symmetric. 

One of the most immediate consequences of the invariance (1) replacing Hermitic- 
ity is that the corresponding Hamiltonian H ^ W admits the existence of complex 
eigenvalues E n ^ E*. In the language of experimental physics this simply means a 
potential instability of the whole system, a phenomenon truly encountered, e.g., in 
relativistic Quantum Mechanics [10] or in certain magnetohydrodynamical systems 



dx 2 



+ V(x) 
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[11]. In the former case the instability may either mean a complete collapse of the sys- 
tem or a mere controlled creation of the particle-antiparticle pairs [12]. An even less 
catastrophic scenario is encountered in experimental magnetohydrodynamics where 
the transitions to the unstable regime prove reversible as a rule [13]. 

The potential theoretical consequences of VT — symmetry seem equally appealing 
at present. One of the reasons is that in eq. (2) the coordinates x themselves may 
be complex and, hence, manifestly unobservable! This means that one gets quite 
close to the relativistic field theory where the corresponding analogue of x (viz., the 
field amplitude ip) is also an auxiliary (i.e., in principle, not necessarily observable) 
quantity [7, 14]. 

In a way inspired by the latter comment the auxiliary role of the complex, un- 
observable coordinate x may equally well be played by any other element x k of any 
other unobservable set V ^ (—00,00). Without getting too deeply involved in this 
theoretical question we decided to parallel the most popular strategy (where V is 
being specified as a suitable complex curve [6, 7]) and to try to treat the (admissible) 
replacement (—00, 00) — > V as a simplification, discretization of the real axis. 

Our present paper summarizes a few of our most interesting results. Firstly, we 
set the scene in section 2 where we list some supplementary motivations and show 
how we intend to construct our illustrative examples. In the subsequent section 3 
both the physical and mathematical background of our considerations is outlined in 
some necessary detail. 

Up to the final brief summary of our results in section 7, the rest of our paper 
becomes more technical and it shows how our simplest discrete models can be solved 
and how their various generalizations can be treated and remain, in a current vague 
sense of the word, exactly solvable. 
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2 PT-symmetry in Quantum Mechanics 



2.1 Non-analytic Schrodinger equations 

The admissibility of a complexification of spectra is a puzzling possibility which in- 
spired qualitative analyses of many concrete models. For many analytic potentials 
V(x) of an immediate physical interest the constructive demonstrations of the sta- 
bility (i.e., of the reality of the spectrum) have been performed by perturbative [15], 
quasi-classical [16] as well as purely numerical [17] methods. 

At an early stage of development the non- analytic solvable models have been 
considered conspicuous. Many of them - typically, potential V(x) ~ \x\ on the whole 
real line, solvable in terms of Airy functions by the matching method - happened 
to exhibit instabilities (i.e., complex energies) at all the non- vanishing values of the 
corresponding coupling constants [8] . This caused a certain delay of return to some of 
the most elementary exactly solvable piecewise-constant VT — symmetric potentials 
in eq. (2). Still, their later studies [18, 19] clarified that their spectrum of bound 
states may remain real in a fairly large domain of their coupling constants (cf. also 
[20, 21]). In ref. [22] it has been proved that our VT— symmetric Schrodinger eq. (2) 
possesses solely the discrete spectrum of real energies E = E n , n = 0, 1, . . . for all 
the complex potentials which are "not too strong" , 

HVloc < ^tt 2 ^ 3.701. (3) 
8 

At the same time, the possible emergence of the instabilities at all the non-zero cou- 
plings has been confirmed for certain piecewise-constant VT — symmetric potentials 
acting on the complex, curved contours [23]. 

2.2 Runge-Kutta option 

The indefiniteness of conclusions of the above observations inspired out present ex- 
tension of the study of the piecewise-constant potentials towards the models which 
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are defined over a mere discrete, equidistant lattice of points 



x = -1, 



x k = Xk-i + h — -1 + kh, 



h 



AT' 



k = l,2,...,N. (4) 



These models possessing VT— symmetry and real spectra will be characterized by 
the variability of the number N + 1 of grid points in (4). We believe that in this 
perspective our understanding of the generic VT— symmetric models may acquire 
a new dimension, returning us to the usual differential-equation framework in the 
N — > oo limit. 

One of the most natural possibilities of the discretization of the differential eq. (2) 
is given by the well known Runge-Kutta recipe [24] 
ip(xk+i) - 2 ip(x k ) + ^(ajfc-i) 



h 2 



Using the standard boundary conditions 



+ V(x k )ip(x k ) = Eip(x k ) . 



(5) 



ip(x ) = i>{x N ) = 

we shall employ here just the piecewise constant, purely imaginary antisymmetric 
potentials of refs. [18, 21], 

+lZ n X E (-£ n , -£ n -l), 

, n=l,2,. + (6) 

i Z n x G (£ n _i,£ n ), 

Their discontinuities lie at the matching points £q — < i\ < . . . < £ q+ i = 1 so that 
at q = we have to solve the discrete Schrodinger bound-state problem 

ip(x k +i) - 2 ip(x k ) + ip(x k -i) 



V(x) 



- isign(x k )Z *jj(x k ) = Ei/j(x k ), 



q = (7) 



etc. It is worth noticing that the first nontrivial N = A version of eq. (7) coincides 
with the Weigert's matrix model 



V 



2 + \iZ -1 





-1 2 -1 
-1 2-\iZ 



\ ( \ 

a 



1 



E 



( \ 

7 



(8) 



studied in ref. [25], in more detail, as "the simplest nontrivial" Runge-Kutta dis- 
cretization of the VT— symmetric square well of ref. [18]. 



3 Physics and mathematics behind the discrete 
Schrodinger equations 



3.1 The definition of the physical metric 

In the light of reviews [4, 9], the PT— symmetric Hamiltonians are not compatible 
with the standard "Dirac's" metric / in Hilbert space (indeed, we have H ^ W). 
A construction of an appropriate (i.e., of a Hamiltonian-dependent) generalization 
7^ / of the metric is necessary for a correct probabilistic interpretation of the 
measurements. This means that our Hamiltonian only must be Hermitian (or, more 
precisely, self-adjoint) with respect to the generalized metric, 



All the physical contents of our models must be derived from the new metric of 
course. For this reason it makes sense to avoid confusion by calling eq. (9) with 
"quasi- Hermiticity" condition [4]. 

One can notice that almost exclusively, standard textbooks train our "quantum 
intuition" on the representations of the Hilbert space where the metric remains triv- 
ial. This means that all the nonstandard PTSQM considerations profit enormously 
from any construction of a sufficiently flexible class of solvable models. In such a 
setting our present finite-dimensional VT— symmetric models may play the role of a 
guide towards the interpretation, especially due to their facilitated connection with 
the continuous N — > oo limit. 

A parallel reason for the attention paid here to the discretized systems lies in the 
well known difficulty of the problem of the efficient construction of 6 ^ I for a given 
H 7^ if* '. Using the Dirac- inspired "brabraket" compact notation of ref. [26] and 
employing the biorthogonal basis composed of the left and right eigenstates of our 
non-Hermitian Hamiltonian H we may recollect (cf., e.g., [27]) that the knowledge 
of the spectral representation of the Hamiltonian, 



ift = QHQ'\ 



I = 0t > 0. 



(9) 



ir = X» 



((n\n) 



n 
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almost immediately leads to the parallel multiparametric formula 

© = EI™»MH, 0n>o (io) 

n 

which specifies the correct metric as a (non-unique, parameter-sequence-dependent) 
solution of the operator eq. (9). Thus, we can employ this formula and construct 
6, more or less comfortably, for the majority of the exactly solvable models of the 
form (2) [28]. This is in contrast with the situation where one must construct the 
states \n)) by some approximative method. This would make formula (10) practically 
useless. Complicated computations would be needed which, typically, search for 6 
in the form of a product where the parity V is multiplied by a charge C (from the 
left, [9]) or by a quasi-parity Q (from the right, [29]). 

3.2 Matching method 

Let us recollect that the N = oo continuous-limit models (6) are exactly solvable 
at any q > and that their practical solution remains feasible and transparent at 
q = [18], at q — 1 [30] and at q = 3 [31] at least. The solvability in closed form also 
characterizes the parallel modifications of these models characterized by the periodic 
boundary conditions [32]. 

It would be highly desirable to extend the available exact solution of the discrete 
Weigert's "minimal" finite-dimensional model (8) with q = and N = 4 to some 
higher integers N and/or q. A key purpose of our present paper is to show that such 
extensions are feasible, indeed, once we adapt the matching method, so efficient in 
its standard differential-equation form, to the needs of difference equations. 

In the next section our constructive considerations will start from one of the most 
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straightforward N > 4 generalizations 

( ii-F -l 

-1 ■■■ 

-1 

-1 



V 



1 



Cir. 



I 



7 



Pn 
V #° ) 



(11) 



Ai-F ••. 

-1 

-1 -i£-F 

of the Weigert's model represented by the N — 1 = 2n + 3 dimensional matrix tran- 
scription of our difference Schrodinger eq. (7) with the re-scaled energy eigenvalues 
F = E h 2 — 2 and with the re-scaled strength £ = Z h 2 of the non-Hermiticity. 

We shall see that such a matrix re-formulation of the matching recipe proves 

solvable in closed form at all n — 0, 1, In the Hermitian context of Quantum 

Chemistry we just arrive at the well known Hiickel's solvable models in the limit 
Z -> [33]. 



4 Finite-dimensional Schrodinger equation (11) 
4.1 Closed solvability 

Recollecting the definition of the classical Tshebyshev polynomials of the second kind 
[34], 

(cf. a collection of their most relevant properties in Appendix A) and using the real 
elements 7 = ^(O), = Reip (xk+i) and bk — Imip (xfc+i) w e are now prepared to 
formulate our first result. 

Theorem 1. 

Whenever the VT— symmetry remains unbroken, closed solutions of eq. (11) are 
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defined by the formulae 



a k = cik + ibk, (3k 
-F + iA 



a k = (a + ib) U k 



a k -ib k = : i 
k = 0,1, 



and 



f-F + iA 

7 = (a + i&) C/„+i = (a - i&) £4+i 



complemented by the selfconsistency constraint 

' -F + iA 



F 7 = -(a + i6) C/ n 



(a - ib) U n 



-F-iC 



-F-iC 



(12) 



(13) 



(14) 



2 J y ' \ 2 
Proof. 

In the first step of our analysis we get the relation (3 k = &l due to the VT— symmetry 
of our eigenvectors. We then identify eq. (11) with the three-term recurrences sat- 
isfied by the Tshebyshev polynomials (cf. Appendix A) and verify formula (12) for 
eigenvectors by showing that it is compatible with the corresponding boundary con- 
ditions (at the smallest subscripts k). As long as all this reduces the full tridiagonal 
(AT — 1) x (N — l)— dimensional matrix eq. (11) to the mere three matching conditions, 

U n-i (^r 1 ) (a + i&) 
U n (a + ib) 

7 



-1 i£-F 
-1 




-1 
-F 
-1 





-1 
-i£-F 



\ 



= 



(15) 



y u u -i -iq - r -l j U n ( F 2 lg ) (a-ib) 

Un-1 (=^) (a - ib) 

it remains for us to show that the first and the third line may be simplified to give 
eq. (13) while the middle line defines the product (14). QED 

Remark. 

Without an assumption of an unbroken VT— symmetry one would have to admit 
that the eigenvalues F are complex. A more explicit illustration of what happens 
has been described, in the continuum limit, in ref. [35]. 

Corollary. 

A nontrivial solution always exists at F — 0. 



Proof. 



In a discussion of consequence of Theorem 1 one has to notice that formula (14) forces 
us to distinguish between two regimes where F = and F ^ 0, respectively. Once we 
set, tentatively, F — 0, it is easy to deduce from eq. (13) that the parameter a must 
vanish for the even n = 0, 2, 4, . . . (and we may normalize 6 = 1) while 6 = and 
a = 1 for the odd n — 1, 3, 5, . . .. Thus, eq. (13) degenerates to the mere definition of 
the last element 7 of the eigenvector and we are left with the single secular eq. (14) 
which acquires the following two alternative forms, 



"guess of the energy" was correct and that F = is always an eigenvalue. QED. 
Remark. 

In spite of the non-Hermiticity of the Hamiltonian, the F = eigenvalue is "robust" 
[30] and remains real at all the real couplings Z e (—00, 00). The spectrum remains 
symmetric with respect to this "middle point" [note that F = corresponds to the 
energy E = E n+2 = 2/h 2 = N 2 /2 = 2(n + 2) 2 ]. The existence of such a "central" 
eigenvalue is not in contradiction with the differential-equation results of ref. [22]. 
Indeed, this level moves up with dimension iV and disappears in infinity in the limit 
iV — * 00. We may conclude that in this sense a "richer" structure is exhibited by the 
spectrum at the finite dimensions. 

4.2 Properties of the "fragile" F ^ solutions 

Whenever F^Owe may treat eq. (13) not only as the condition of vanishing of the 
imaginary part of 7, 



U n {\li)-U n {\\i)=0, 

u n (\i£) + u n ( r \it) = o, 



n = 2m + I. 



n = 2m, 



(16) 



These conditions are satisfied identically at any m — 0, 1 



Thus, our tentative 




(17) 
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but also as an explicit definition of the non-vanishing left-hand-side quantity F7 in 
eq. (14). Its insertion simplifies the latter relation, 

T n+ i (-^) (« + ib ) = (-^) (° " ^ ( 18 ) 

where T k (z) denotes the k— th Tshebyshev polynomial of the first kind. 

One of the latter two relations defines the normalization vector (a, b) = (a , 60) 
while their ratio gives 

^ (=*±*) u, +1 (z^ii) + r „ +1 (z£^) ^ (z£±M) . 0. m 

This represents our final secular equation which defines, in an implicit manner, the 
energies F as functions of the couplings £. 

An efficient numerical treatment of the latter eigenvalue problem may be based 
on the re-parametrization 

—F + i£ 

= cos ip, Ref = a, lm<p = (3 (20) 

i.e., 

- F = — cos a cosh (3, - £ = — sin a sinh (3 (21) 
and, in the opposite direction, 



1 



cosa = ^-;F, sinh/3 = — L \ F 2 + £ 2 - 4 + J (F 2 + £ 2 - 4) 2 + 16£ 2 . 

2 cosh /3 2^ v v 

This change of variables transforms eq. (19) into the trigonometric secular equation 

Re sin[(n + l)y9] cos[(n + l)y?*] = Q 
sin (p 

Its roots must be determined, in general, numerically. 
Lemma. 

In the domain with negative (3 < the roots of eq. (22) have a G (0, 7r/2) at the 

negative F < and a G (7r/2, 7r) at the positive F > 0. 

Proof. 

The constant value of the coupling £ > is mapped upon a downwards-oriented half- 
oval in the a — (3 plane. Its top lies at a = 7r/4 while its two asymptotes a = and 
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a — 7r/2 are reached in the limit (3 — > — oo. In such a representation the "robust" and 
£— independent energy level F = lies on the top of the half-oval while its decreasing 
and increasing neighbors are found displaced to the left and right, respectively, along 
the half-oval downwards. QED. 

Remark. 

At the first few lowest dimensions N — 1 = 2n + 3 the roots of eq. (22) may be 
written in closed form, 

F = 0, F ± = ±v/2^F, n = 0, 
F = 0, F ±!± = ±^2 - e ± V^W, n=l 

etc. This enables us to determine the respective critical values Z = Z^ crit ^(N), i.e., 
the exact Z crit (A) = 4 \[2 & 5.66 at n = and Z crit (6) = 9/2 = 4.50 at n — 1 followed 
by the numerically calculated values Z crit (8) = 4.463 at n — 2, Z crit (10) = 4.461 at 
n = 3 and Z crit (12) = 4.463 at n = 4 etc. These results do not contradict the 
expected n — > oo limit as derived in ref. [35], Z crit (oo) = 4.475 (cf. Figure 1). 



5 The role of the integers TV and q 
5.1 Closed solvability at the odd N = 2n + 3 

In comparison with ref. [25] a "one-step easier" discretization of eq. (2) emerges at 
q = and at the odd N = 3, with the two energy roots. A non-equivalent discrete 
alternative to eq. (11) could then read, in the same notation, 



if - F -1 
-1 '<-F ••. 

-1 


-1 
i£-F 


-1 




-1 


-i£-F ••. 






-1 






-1 -if-F 



a. 



a,. 



0. (23) 
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This leads to the following alternative result. 
Theorem 2. 

Whenever the VT — symmetry remains unbroken, closed solutions of eq. (23) are 
defined by formulae (12) accompanied by the alternative matching condition 



(a + ib) U n+1 = (a - ib) U n (-!-— * 



(24) 



(25) 



and by the odd— N counterpart of eq. (19), 
Proof. 

Definition (12) of the eigenvectors remains unchanged but the matching condition 
generated by the subproblem 



V 



-1 i£-F -1 
-1 -i£-F -1 



C/n-i (a + ib) 

U n (a + ib) 

Un l=^) (a-ib) 

U n -! (=^) (a - ib) 



(26) 



is just one, a n+ \ = a* n . The ratio between this equation and its Hermitian conjugate 
eliminates all the normalization ambiguities and leads to the secular equation (25). 
QED. 

Remark. 

Secular equation (25) is an implicit definition of the N = 2n + 3 energy levels 
F = F(£). It possesses the compact analytic solutions at the first few values of n of 
course, 



F± = ±VT=e, n = 0, 
F ±i± = ±\ v / 6-4e 2 ±2 v ^T6F, n = 1. 
The respective elementary expressions for the critical constants 

Z cri t(3) = 1 = 2.25, n = 0, 



Zcrit(5) = ^fl » 3.49, 



n = 1 
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are followed by the complex Cardano representation of the real Z crit (7) = 3.946 at 
n = 2. At n > 2 one switches to a purely numerical algorithm giving Z crit (9) = 4.148 
at n = 3 etc. In comparison with the parallel results derived at even N in section 4 
above we notice that the numerical convergence towards n = oo is perceivably slower 
and from the opposite side at the odd iV (cf. Figure 1 once more). 

5.2 Models with more matching points 

New features of the spectrum emerge when we choose q > in eq. (6). In the light 
of the differential-equation results [30, 31] one may expect the existence of the so 
called 'robust' levels at q > 1 and iV < oo. By definition they remain real at the 
arbitrarily large coupling constants Z k . A priori, the position of these levels may be 
expected to be controlled by both the parameters Zk and £k- 

Obviously, the interpretation of the discretization becomes facilitated when we 
choose the discontinuities at the simple rational numbers £k- It is interesting to note, 
marginally, that their choice simplifies even the trigonometric secular determinants 
in many differential-equation models. Nevertheless, in our present context it is much 
more important that the choice of the simplest £ k at q > 1 enables us to continue 
working with the not too large matrix dimensions N — 1. 

5.2.1 I = 1/2 and iV = 6 

For the q = 1 interaction model (6) with l\ = 1/2, 





= < 







for 



(27) 
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we may start from its discretization (5) with the most elementary choice of iV = 6, 



i£-F 


-1 




-1 


—F -1 
-1 —F -1 
-1 —F 


-1 




-1 


-i£ - F 



7o 
7 
To 



0. 



(28) 



This equation determines the following five eigenvalues, 

1 



F = 0, F ±)± = ±-y 8 - 2 e ± 2 ^4 + £ 4 - 

Besides the constant and £— independent level Fq = one can show that the two 
"external" levels 



±v / 3 



i - - + — e 4 + — + o (e) 



12 ' 288 3456 
remain robustly real for all the values of £ and never complexify, 

31 



±^2 



1 9 1 4 

1 + -V V 

A y 32 y 128 



2/ 6 + O (y 8 ) 



They are complemented by the fragile pair of levels F±- which coincide with F = 



at the critical value of ^ crit = y3/2 pa 1.2247 and complexify beyond this "excep- 
tional" point. 



5.2.2 



1/2 and TV = 8 and more 



The N = 8 calculations lead to the observation that the two outmost energy levels 
as well as the "middle" F = level remain robust and real at all the real £. At 
the critical value of ^ crit = 0.845479352 we observe the confluence and subsequent 
complexification of both the remaining (i.e., positive and negative) energy pairs at 
the respective two numerical values of F crit = ±1.05167218. With the further growth 
of £ the quantity Re F 2 becomes negative beyond £ 2ero = 3.2222152046, with the 
purely imaginary crossing-point F zero = ±2.146638195 i. 
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At N = 10 the secular algebraic equation for the energy-level characteristics 
y = F 2 has the following four roots in the Hermitian limit with £ = 0, 

5 1 r- 3 1 r- 

y a ,± = - ± - V5, y b: ± = - ± - V5. 

Numerically one finds that the real quadruplet 

0.5173571919, 1.810807242, 1.810964520, 3.356678112 

of these values gives the single roots F = ±0.7192754632 and F = ±1.832123935 
and the double roots F = ±1.345662380 at £ subcri t = 0.50209209. In parallel, the 
complex quadruplet 

0.5173571921, 1.810885878 ± 0.00005268889724 i, 3.356678117 

emerges numerically at the upper estimate of the critical coupling strength £ snperC rit = 
0.502092091. 

Similar observations have been made also at iV = 12 and iV = 14 and all of them 
parallel and reproduce the qualitative features as observed in ref. [30] in the iV — > oo 
continuum limit. 



6 The role of the shift-parameter t 



6.1 Suppressing the non-Hermiticity using I = 5/8 



In order to illustrate the emergence of a pair of the robust levels let us re-consider 
the q — 1 interaction model (6) and employ a larger shift t\ — 5/8. In the resulting 
potential 



V(x) 



+iZ 


-iZ 



for x G < 



1,- 

5 5 



8) ' 



(29) 



1 (1.1) 
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the non-Hermiticity is weakened at a fixed coupling Z. We get the matrix analogue 
of problem (28) from the discretized eq. (5) with N = 8, 



i£-F 



V 



-1 



— F -1 

-1 -F -1 

-1 — F -1 



-1 — F -1 
-1 -F 



«0 



-i^-F ; 



7i 
7o 
7 
7o 

V a *o J 



(30) 



An inspection of the explicit form of the secular determinant of this equation, 



V 



-F 6 - F 4 (e 2 - 6) + F 2 (4£ 2 - 10) - 3£ 2 + 4 



(31) 



confirms that all of its seven roots remain real in the Hermitian £ — > limit, 

F = 0, F ±i0 = ±v / 2, F ±)± = ±V / 2±v / 2, £^0. (32) 

Moreover, in a way illustrated in Figure 2, only two roots F± _ merge and complexify 
beyond £ cr ^ ~ 1.15470. In the other words, as many as five of the roots remain real 
for Z > 1, 

F) = 0, F ± ,o = ±l, F ±i+ = ±v / 3, £^cx). 

This observation parallels the similar results obtained in the continuous limit. All 
this also agrees with the a priori expectations, reflecting the weakened influence of 
the non-Hermitian part of the potential at the comparatively large shift £ = 5/8. 



6.2 Strengthening the non-Hermiticity using £ = 3/8 

The number of the robust levels may vary with the width I of suppression of the 
imaginary interaction. Using a weaker suppression t\ = 3/8 we may modify eq. (29) 
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accordingly and we get the precise N = 8 analogue 



V 



—1 


-1 


— 1 












-1 


—F 
-1 


-1 
—F 
-1 


-1 
—F 


-1 












-1 


-i£-F 
-1 


-1 
-i£ - F 



«1 



7o 
7 
7o 



(33) 



of problem (30). It is easy to evaluate the related secular determinant 



V 



_ F 6 _ F 4 ^2 _ ^ + F 2 ^4 + 4 ^2 _ + 2 ^4 + ^2 + 4 



(34) 



which gives the same roots (32) in the weak-non-Hermiticity-limit £ — > 0. A change 
may be expected to emerge in the more non-Hermitian regime. Quantitatively this 
change is being well illustrated by Figure 3 where we see that merely two non- 
vanishing roots (viz., F± :+ — >• ±y/2) remain robust and real at large Z — > oo. We see 
that two complexifications take place at a certain finite coupling Z, i.e., at a unique 
critical point £ cr j t ~ 0.5875691807. In contrast to the previous case we witness two 
distinct mergers of the eigenvalues F± t0 with F±- at the two different exceptional 
points Fl xc - w ±1.140716421. 



7 Summary 

A tentative weakening of the standard Hermiticity H = W to the mere V— pseudo- 
Hermiticity of the Hamiltonian, 

H^=VHV- X (35) 

employs only too often the ordinary differential Hamiltonian H and the parity op- 
erator V [7]. Here we intended to emphasize that it is by far not the only possible 
scenario. For this reason we replaced eq. (2) by its discrete counterpart (5) and 
explained that its use may bring a few important advantages. 
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We feel satisfied by the fact that the purely formal application of the PTSQM 
recipe survives quite easily our present narrowing of one's attention to the very spe- 
cific and solvable square-well models. We have shown that the use of various specific 
discrete versions of the Hamiltonians leads to very relevant formal simplifications 
emphasized throughout the text. 

Within the dynamical region where our models possess the real spectrum, very 
close parallels were shown to emerge between the bound states for the standard 
continuous coordinate x [18, 30, 31] and for its present less standard discrete coun- 
terpart (s). 

We have seen that in the discrete case it is particularly easy to complement the 
validity of the auxiliary, technical pseudo-Hermiticity condition (35) by the construc- 
tive solution of the vital and essential quasi-Hermiticity requirement (9) which is of 
key relevance in the context of physics. 
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Figure captions 

Figure 1. Numerical values and convergence of the critical 
couplings Z {cnt) = Z (cnt) (N). 

Figure 2. The £— dependence of the seven roots of the sample 
secular determinant (31) 

Figure 3. The £— dependence of the seven roots of the sample 
secular determinant (34) 
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Appendix A: Tshebyshev polynomials of scalars and matrices 



A "pedestrian's" way to the simultaneous work with the Tshebyshev's polynomials 
Q n (x) of the first kind [replace Q n (x) — > T n {x)] and of the second kind [replace 
Qn(x) U n (x)] leads through the "universal" three-term recurrence relation 

Q n+1 (x)-2xQ n {x) + Q n . 1 (x) = (36) 

complemented by the specific initial values, 

U (x)=0, U 1 (x) = VT^, ... 

T (x) = 1, 7i(ar) = x, .... 

Many useful identities valid for these polynomials are available in the literature. 

For complex arguments, an interesting situation emerges during a real-matrix 
re-arrangement of eq. (11). We may split eq. (11) in its real and imaginary parts and 
'glue' them together in a pentadiagonal eigenvalue problem 

a 
bo 



-F -Z 

z -f 


-1 
-1 








-1 
-1 


-F -Z 

Z ~F 














-1 
-1 








-1 
-1 


-F -Z 
Z -F 


-1 









-2 


—F 



V 

It can be written in a partitioned form, 
'X -1 

-l x ••. 
■•. ■•. -1 

-1 X 



X 7 / 



= 0. 



(37) 



Ci 



2d T —F 
24 



= 0, 



(38) 



V 7 J 



with an 'odd' anomalous row and column containing an auxiliary two-dimensional 
vector [d T = (1,0)]. 

In the new context all our wave functions are again proportional to the Tsheby- 
shev polynomials, 



1. 



Ck = U k [ -X c , 



X = 



k = 0,l,...,n + l. 



(39) 



-F -i 

V i -f ) 

The existence of the explicit solutions (39) [where polynomials Uk depend on a two- 
by-two matrix argument X] reduces eq. (38) to the constraint imposed upon the 
two-dimensional real vector c n+ i. Of course, this condition is equivalent to the 
complex matching as mentioned above. 
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